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ABSTRACT 

We use deep wide-field F-band imaging obtained with the Wide Field Camera at 
the prime focus of the Issac Newton Telescope to study the spatial and luminosity 
distribution of galaxies in three low redshift (0.04 < z < 0.2) clusters: Abell 119, 
Abell 2443 and Abell 2218. The absolute magnitude limits probed in these clusters 
are My — 51og/i .7 = —13.3, —15.4 and — 16.7mag, respectively. The galaxy pop- 
ulation, at all luminosities, along the line-of-sight to the clusters can be described 
by the linear combination of a King profile and a constant surface density of field 
galaxies. We find that, for these three clusters, the core radius is invariant with in- 
trinsic luminosity of the cluster population to the above limits and thus there is no 
evidence for luminosity segregation in these clusters. The exception is the brightest 
galaxies in A2218 which exhibit a more compact spatial distribution. We find the to- 
tal projected luminosity dis tribution (within l/i^Mpc of the cluster centre) can be 
well represented by a single iSchechterl l)1976f) function with moderately flat faint -end 
slopes: a = -l.22±^% (A119), a = — l.ll±g;Jg (A2443) and a = -1.14±g;[$ (A2218). 
We perform a geometric deprojection of the cluster galaxy population and confirm 
that no 'statistically significant' evidence of a change in the shape of the luminosity 
distribution with cluster-centric radius exists. Again, the exception being A2218 which 
exhibits a core region with a flatter faint-end slope. 
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1 INTRODUCTION 



Galaxies of different types in clusters are known to have 
dif ferent projecte d spatial distributions. This was realized 
bv lOemieil jl974f) . who showed that less luminous galaxies 
have a more extended profile th an t he more massive ellip- 
ticals. iMelnick fc SargenrJ (Il977l) and iDresslerl il98Ct> iden- 
tified what is now known as the 'morphology-density' re- 
lation, where the relative fractions of elliptical, lenticular 
(SO) and spiral galaxies depend on the surface density, while 
IWhitmore et alJ <ll993l) argued that these trends are better 
correlated with cluster-centric radius. In a recent compre- 
hensive stu dy of an ensemble cl uster built from 59 nearby 
rich clusters lBiviano et al.l (120021) . demonstrated clear segre- 
gation between ellipticals, early and late-type s pirals. This is 
also seen in a single HST mosaic of Abell 868 bv lDriver et alJ 



(2003) which concludes that cluster cores are devoid of, or 
at least depleted, in late-type systems. 

However, as well as morphological segregation, ev- 
idence is also emerg ing for luminosity segregation. 
iRood fc Turnrosel lll96St) argued that dwarf s were less con- 
centra ted than giants in the Coma cluster; ICaoelato et alJ 
ll98l|) detected a mass-density relation in Abell 196; 
lYepes et al.l 11199 If) studied luminosity segregation in a num- 
ber of clusters and found that the degree of segregation 
correlates with the dynamica l state of the cluster. The 
study of iFerguson fc Sandaeel il989T) in Virgo and For- 
nax demonstrated that dwarf ellipticals were highly con- 
centrated leading to a division of the dwarf population 
into distinct strongly clustered nucleated dwarf ellipticals 
and a distributed population consisting of non-nucleated 
dwa rf ellipticals and dwarf irregulars. In th e Com a clus- 
ter, lLobo et all il997t) and iKashikawa et alJ ll 19981) found 
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evidence for strong luminosity segregation, with the giants 
being clumped in two substructures whi le the dw a rfs tra ced 
a more diffuse and regular distribution. lAndreonl 1)20021) ar- 
gued that some form of mass segregation is also at work in 
the z = 0.31 cluster AC 118 (also known as Abell 2744). 
The giant ellipti cals and len ticulars may also be kinemati- 
cally segregated (ISteinlll997l) . suggesting that these objects 
are the original kernel of the clusters while spir als and dwarfs 
are co mparatively late arrivals. Conversely, iBiviano et alJ 
( 2002) find that the only evidence for luminosity segrega- 
tion is for ellipticals outside of substructures in their ensem- 
ble cluster. 

ISmith et all Jl997t) and lDriver et alJ <1998h found that 
the dwarf-to-giant ratio shows a trend with density and, be- 
cause of the approximately spherical shape of clusters, with 
radius. This lead directly to the idea of a dwarf-density re- 
lation |P hillipps et ftl-||l99 ^- Together w ith previous work 
llLobo et^dTEj99«lKasMkawa et al]ll99s|) this may suggest 
that dwarfs are especially affected by the cluster environ- 
ment, as one would expect for such fragile objects. 

IPracv et al.l i2004h have recently investigated luminos- 
ity segregation using a wide, deep mosaic of HST images 
of Abell 2218 and found evidence that dwarf galaxies avoid 
the central regions of this cluster and trace a more spa- 
tially extended distributi on. A similar result wa s found for 
the NGC 5044 group by iMathews et all J2004h and for a 
sample of loose groups bv lGrrardi^^dl J2003ft . If the segre- 
gation for dwarfs is real, it may originate from initial condi- 
tions, where low luminosity galaxies are only now in-falling 
into clusters (e.g., Croton et al. 2005), or it may be due to 
processes internal to clusters, such as tidal disruption and 
galaxy harassment. For these reasons, it is important to in- 
vestigate the existence of luminosity segregation for dwarfs 
in a broader range of objects and to study its correlation 
with cluster properties. This is now feasible by wide-field 
imaging of relatively nearby clusters with panoramic mosaic 
cameras on 2m telescopes, and we present here the results 
of such a study for three clusters observed from the Isaac 
Newton 2.5m Telescope with the Wide Field Camera. 

In this paper we use relatively-deep wide-field imaging 
of three galaxy clusters in the redshift range 0.04 < z < 0.2 
to measure the galaxy surface density in intervals of lu- 
minosity, recover the overall luminosity distributions and 
the deprojected luminosity distributions to explore the spa- 
tial segregation of the galaxy population with luminosity. 
The absolute magnitude limits probed in these clusters are 
My — 51og/io.7 = —13.3, —15.4 and — 16.7 mag. In Section 
2 we introduce the data and describe its reduction and the 
detection analysis and classification strategy. In section 3 we 
fit a King profile plus constant offset in intervals of luminos- 
ity to simultaneously determine both the cluster population 
profile and the non-cluster foreground/background level. In 
Section 4 we derive the projected luminosity distributions 
in these clusters and in Section 5 we perform a geometric 
deprojection to recover the 'true' luminosity distribution of 
the galaxy population. We summarise our findings in Sec- 
tion 6. Throughout we adopt a Am = 0.3, = 0.7 and 
Ho = 70kms _1 Mpc -1 cosmology. 



2 THE DATA 

2.1 The Observations 

The sample consists of three Abell clusters: A119 (Rich- 
ness=l, BM=II-III, z=0.044), A2443 (Richness=2, BM=II, 
z=0.108) and A2218 (Richness=4, BM=II, z=0.181). The 
observations were obtained on the nights of 2 nd & 3 rd 
September 2000 using the Wide Field Camera (WFC) 
mounted at the prime focus of the Isaac Newton Telescope 
(INT). The WFC consists of a mosaic of four 2048 x 4096 
thinned EEV CCDs with a plate scale of 0.333 arcsec/pixel. 
The total sky coverage is 0.287 deg 2 per pointing. The imag- 
ing of each cluster consists of four partially overlapping 
pointings with the WFC, the exception being the higher 
redshift cluster A2218 which is a mosaic of just 2 pointings. 
Each cluster was imaged through the V filter with an expo- 
sure time of 1200s. A summary of the observations is given 
in Table □ 



2.2 Data reduction 

The data reduction was performed by the Cambridge Astro- 
nomical Survey Unit (CAS U) and full detai l s of th e pipeline 
procedure can be found in llrwin fc Lewisl i200ll) . In sum- 
mary, the data are first bias subtracted and trimmed. Bad 
pixels and columns are interpolated over using data from 
neighbouring regions. All four chips are then corrected for 
non-linear behaviour in the two analogue-to-digital convert- 
ers. The data are then flat fielded using master sky flats and 
a gain correction is applied so that all the CCDs have the 
same zero-point. Finally, an astrometric solution is derived 
by matching bright stars in the field-of-view of each chip to 
the Guide Star Catalog. 



2.3 Photometric calibration 

Four lLandold lll992fi standard star fields (SA92, SA95, 
SA110 and SA113) were observed at various air-masses 
throughout each night. For each observation of a standard 
star a zero-point was calculated: 

ZP star = m + 2.5 log I (1) 

where m is th e magnitude of the standard star given by 
lLandold dl992l) . / is the star counts in ADUs as measured via 
SExtractor fixed aperture magnitude (7" radius) and t is 
the exposure time (10s). These zero-points were fitted with a 
dou ble-linear funct ion in airmass and (B — V) colour (taken 
from lLandoldll992T) to derive the extinction coefficient (a a ir) 
and an above-atmosphere zero-point ao. 

ZPstar = ao + a air (sec Z) + a co i (B — V) (2) 

We then assigned to each field a final zero-point ZP = 
ao + a a irsec Z. Fig. shows the data used for the ^-band 
calibration on the night of 3 rd September 2000, the top panel 
shows the colour-corrected zero-point versus airmass and the 
bottom panel shows the airmass-corrected zero-points versus 
(B — V) colour. The fit (equation to the data is repre- 
sented by the solid line and the RMS of the data (top panel) 
is ±0.06 mag. 
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Table 1. Data characteristics. 
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Figure 1. V-band standard star calibration for 3 Septem- 
ber 2000. Colour-corrected zero-points versus airmass (top) and 
airmass-corrected zero-points versus colour (bottom) 



subset of detected objects in this plane, for each cluster, is 
illustrated in Fig. El The central surface brightness was cal- 
culated in a circular aperture with an area equal to that of 
the detection criterion (i.e., 8 pixels). 

2.6.1 Cosmic ray rejection 

In Fig. El a group of objects with a high central surface 
brightness (at a given magnitude) is clearly discernible (up- 
per right). These objects, which have surface brightnesses 
higher than that of stars, are cosmic rays. We therefore de- 
fine a region in the jj, -V plane such that: 

fio < aV + b (3) 

and classify all objects in this region as cosmic rays. Equa- 
tion © is shown as the dot-dashed line in Fig. [5] The slope 
(a) and intercept (b), were chosen separately for each clus- 
ter, to best match the data. The objects classified as cosmic 
rays are shown as blue points in Fig. |5] 



2.4 Object detection and photometry 

Objects were detected automatically u sing the SExtrac- 
TOR package feertin fc Arnoutslll996f) . SEXTRACTOR de- 
tects objects as groups of connected pixels which exceed a 
certain threshold. A criterion of 8 connected pixels above an 
isophotal threshold of fio = 26.0 magarcsec -2 was used. The 
best_mag parameter was used to measure total galaxy mag- 
nitude (hereafter denoted V) ; this corresponds to a iKronl 
( 1980) extraction aperture, except for crowed regions where 
an extrapolated isophotal magnitude is used. The Kron ex- 
traction aperture is set to 2.5 times Rk where Rk is the first 
moment of the image distribution. In regions where point- 
ings overlap, the duplicate objects were removed from the 
catalogs. 

2.5 Exclusion regions 

After the initial object detection each field was visually in- 
spected and any CCD defects, very bright stars, diffraction 
spikes and satellite trails causing spurious detections were 
identified. Circular and rectangular regions enclosing these 
areas were defined and excluded from the catalogs. In ad- 
dition, the areas within 20 pixels of the CCD edge and the 
vignetted corner of CCD3 were also excluded. 

2.6 Object classification 

We used the position of objects in the central surface 
brightness-magnitude (n -V) plane to classify objects as 
either galaxies, stars or cosmic rays. The distribution of a 



2.6.2 Star-galaxy separation 

The process of star-galaxy separation begins by identify- 
ing saturated stars in the catalogs. Stars brighter than 

V « 17.2 mag are saturated, these are clearly identifiable 
in the fi -V plane as a horizontal locus of points with 
fi < 18.0 mag arcsec -2 . We classify objects with: 

Ho < 18.0 and V < 17.2 (4) 

as flooded stars (horizontal dashed line in Fig.|2J. The stellar 
locus can be seen in Fig. [5] as a diagonal locus of points 
(red) with a higher surface brightness than the overall galaxy 
population (black points), and extending from V w 17.2 to 

V ~ 21 mag. We therefore define a line in the \i -V plane: 

Ho = aV + b' and V > 17.2 and V < 21.0 (5) 

to pass between these populations, and we use it as a divider 
to separate stars and galaxies (see diagonal dashed line in 
Fig.^J. All objects with V < 21 mag which are classified as 
stars are shown in red in Fig. 

For objects fainter than V ~ 21 mag star-galaxy sepa- 
ration becomes problematic. At these magnitudes the stellar 
locus merges with that of the overall galaxy population and 
the two can no longer be distinguished. In order to per- 
form star-galaxy classification faintward of V = 21 mag w e 
use a similar method to that outlined in lLiske et alJ (2003). 
Since the (logarithmic) slope of the st ar counts should re- 
main roughly constant to V ~ 24 mag feiimmel fc Waenerl 
l200ll) . we can use the star counts derived from the bright 
objects in the catalog and extrapolate them to derive the 
expected number of star counts at fainter magnitudes. We 
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Figure 2. Distribution of detected objects in the central surface 
brightness— magnitude plane. Blue points are objects classified as 
cosmic rays, red points are objects classified as stars directly from 
their positions in this plane, green points are objects classified as 
stars by extrapolation from the brighter star counts and black 
points are the objects classified as galaxies. Top panel: A2218. 
Middle panel: A2443. Lower panel: A119. Only one in six objects, 
randomly selected, are displayed. 



then classified objects as stars based on their position in the 
Ho~V plane (objects with the lowest value of fi — mV) until 
we had obtained the predicted number of stars. Although 
this method will result in some individual objects having 
incorrect classifications, the overall statistical proprieties of 
the catalogs will be correct. The objects that have been clas- 
sified as stars in this way are shown in green in Fig. H 



2.6.3 Completeness 

At this point every object in the catalogs has been classi- 
fied as either a star, galaxy or cosmic ray. We are limited to 
galaxies with a central surface brightness (over an area of 8 
pixels) of fi < 26 mag arcsec" 2 and from Fig. [5] we see that 
objects with central surface brightnesses close to this limit 
only occur in significant numbers at V > 23.5 mag, indicat- 
ing t he beginning of detection incompleteness iGarilli et alJ 
1999). We therefore, define an apparent magnitude limit of 
V = 23.5 mag. 



3 RADIAL PROFILES 

The wide field-of-view provided by the WFC mosaics enable 
us to survey the radial distribution of galaxies beyond the 
domain of the cluster, well out into the surrounding field. 
To explore this, we plot for each cluster the galaxy-surface 
density against cluster-centric radius. This is achieved by 
deriving the galaxy counts in concentric annuli centred on 
the brightest cluster galaxy. The surface density was calcu- 
lated taking into account the area of each annulus which en- 
compasses the unmasked field-of-view of the available CCD 
area (i.e., corrected for any exclusion regions which inter- 
sect the an nulus). We also c orrected for the 'diminishing 
area effect' iDriver et alJ l9981 whereby the area over which 
faint objects can be detected is reduced by the presence of 
brighter objects. To calculate this effect we used the SEx- 
TRACTOR parameter ISOAREA - which returns the total 
area assigned to an object by SExtractor - to calculate 
the amount of area occupied by brighter objects. 

The SExtractor best magnitudes were corrected 
for ga lactic extinction using the maps of ISchleeel et alJ 
(1998). The observed radial galaxy surface density profiles 
for A2218, A2443 and A119 are shown in Fig. EH Fig. Hand 
Fig. H respectively, as open squares. The error bars on the 
points are those expected from purely Poisson statistics. 



3.1 Cluster profiles 

The radial profiles in Figs. HE1 represent the superposition 
of the cluster surface density profiles and a constant 'field' 
galaxy surface density. We elect to descr i be the cluster sur- 
face density in functional form by a lKind il962T) profile, plus 
a constant surface density of 'field' galaxies, thus: 



cr(r) 



(To 



1+ T- 



+ N f 



(6) 



where o~(r) represents the radially (r) dependent number- 
counts along the line-of-sight. The first term in equation 
© represents the projected distribution of cluster galaxies 
and the second term the superimposed 'field' population. In 
Pigs. |2JEl we have binned the galaxy number counts in radial 
annuli. We therefore need to integrate our fitting function 
(equation over the annuli to obtain the average surface 
density, which gives 



2nrcr(r)dr 



• [jV/r 2 -r- ff0 ^log(14-(r/r c ) 2 )] 



(7) 
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Figure 3. The surface density of galaxies in A2218 as a function of cluster-centric radius for a series of magnitude intervals. The solid 
curves are the best fitting 'King + constant' profile given by equation 0. The (apparent) magnitude intervals and their corresponding 
absolute magnitude intervals (for the cluster redshift) as well as the best fitting values of the fit parameters are tabulated in each panel. 
The dashed lines show the value of the Nf parameter from equation 171 . 



where r m i„ and r ma x are the inner and outer radial bound- 
aries of the bin, respectively. The fitted profiles to the mea- 
sured galaxy surface density profiles are shown as the solid 
lines in Fig.|HJ|5]in 1 mag intervals of V. Overall, the surface 
density of galaxies are well described by equation ©. 

In Fig. HJwe show the King profile core radii [given by 
the parameter r c in equation (J7J] as a function of absolute 
1^-band magnitude, the errors on the points are calculated 
directly from the covariance matrix returned from \ 2 min- 
imisation. In the case of the clusters A2443 (middle panel) 
and A119 (lower panel) we find that the core radius is es- 
sentially independent of magnitude, with both clusters hav- 
ing core radii of r c ~ 0.2-0.5 Mpc at all luminosities. The 
best fitting slopes are given by —0.001 ± 0.033 Mpc mag -1 
and 0.058 ± 0.040 Mpc mag -1 , respectively - both consis- 
tent with zero at ~ la. In A2218, we find that the bright- 
est galaxies have a smaller core radius then their fainter 
counterparts, with marginal evidence that this trend - in- 
creasing core radius with decreasing luminosity - continues 
for the intermediate population. The best fitting slope is 
given by 0.086 ± 0.030 Mpc mag -1 . However, if the bright- 
est points (My < —21 mag) are removed, the slope becomes 
0.051 ± 0.054 Mpc mag -1 which is consistent with zero at 
less than the la lev el. This trend is generally consistent with 
IPracv et alJ i2004h who found that the spatial distribution 
of galaxies in A2218 is more extended for the lower lumi- 
nosity populations. Unfortunately the data do not extend 
to their dwarf and ultra-dwarf regimes. 
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Figure 6. Core radii [r c from equation 0] versus absolute V- 
band magnitude. Top panel: A2218. Middle panel: A2443. Lower 
panel: A119. The best fitting slopes are displayed as solid lines. 
The 'Giant', 'Dwa rf and 'Ultra-dwarf regimes, as defined by 
are delineated by the dashed lines. These def- 
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initions have been adjusted due to differences in the filters and 
the assumed cosmology. 



3.2 Reference field counts 

In order to study the cluster galaxy population we need 
to remove the contribution to the counts along the line-of- 
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Figure 4. Same as Fig. [3 except for A2443 
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sight from foreground and background 'field' galaxies. To 
do this we use the standard technique of statistical field 
subtraction. However rather than use data obtained from off- 
cluster pointings we are now able to use the counts derived 
from the radial fits (i.e., Nf in equation . 

The field galaxy number counts derived in this way are 
shown in the top panel of Fig. [7] where we have plotted 
the counts from the A22f8 field (filled squares), the A2443 
field (filled circles) and the A119 field (stars), separately. 
For comparison we also plot the galaxy nu mber counts 
from t he Millennium Galaxy Catalog (MGC) of lLiske et alJ 
(2003), which we have converted from the B-band to the V- 
band using the mean field galaxy colour (B — V) = 0.94 
iNorberg et alJ 120021) . We note that this colour is only 
strictly valid for V < 18 mag; however, the number counts 
agree quite well over the entire range of luminosities. The 
'field' counts derived for the low redshift cluster A119 gen- 
erally have larger error bars and much larger scatter than 
those derived for the other two clusters. This is expected 
since the radial profiles for this cluster extend to a radius of 
only ~ 2Mpc. The scatter in the points can be seen more 
clearly in the lower panel of Fig. |7|where we show the galaxy 
counts relative to a N oc io 0- relation. 

In order to give a smooth representation of the 'field' 
counts we fitted a quadratic function to them. Since 
'field' galaxy counts can vary significantly between fields, 
a quadratic fit was performed separately for each set of ref- 
erence 'field' counts. These fits are shown in the lower panel 
of Fig.|7|as the solid, dashed and dot-dashed lines for A2218, 
A2443 and A119, respectively. For illustrative purposes a si- 
multaneous fit to all three data sets is shown as the solid line 
in the top panel of Fig. [7| We also used the MGC counts at 
V < 18 mag to provide a bright end 'anchor' to the fit - 
we artificially introduced errors of 10% on these points so 
as not to allow the MGC points to influence the fit at the 
critical faint end. The quadratic fits yield: 



dN 
dV 

a 
b 

c 



_ ^Qa + bV+cV 

= -7.83486, -7.64871, -10.87945 
= 0.6404656, 0.6162334, 0.9579601 
= -0.0052439445, -0.0044137863, 



(8) 



-0.0130500505 



(for A2218, A2443, A119, respectively). We use these ana- 
lytical representations of the 'field' counts, adjacent to each 
cluster, for statistical field subtraction throughout. We note 
that this method will also correct, in a statistical sense, for 
any incompleteness in the cosmic ray rejection. 

Since we will use the analytic representation of the 
'field' counts given in equation © to remove contamination 
by 'field' galaxies superimposed on the cluster, we require an 
estimate of the uncertainty in this analytic representation. 
To do this we use the formal variances and covariances for 
each of the fit parameters (a, b and c) returned from the x 2 
minimization. The error (e) on the 'field' counts is given by: 



(5l + V 2 Si + V A 5 2 C 



+2V5 ab + 2V 2 5 ac + 2V 3 S b 



(9) 




Figure 7. Top: 'field' galaxy counts derived from equation J7| 
see text for details); filled squares are the reference 'field' counts 
from the A2218 field; filled circles for the A2443 field and stars 
for the A119 field. The open triangles are the 'field' galaxy counts 
from the MGC. The solid line is a quadratic fit to the counts from 
all three clusters. Bottom: 'field' galaxy counts shown relative to 
a 10 0A5V relation. The quadratic fits to the counts are shown 
as the solid line (A2218), dashed line (A2443) and dot-dash line 
(A119). 



8 2 , 8 2 and 8 2 are the variances and 5 a b, S ac and &bc are the 
covariances on the fitted parameters a, b and c. 



where is the number of galaxies per magnitude interval 
[given by equation JHj] centred on apparent magnitude V . 



3.3 Lensing 

A possible concern regarding the background counts is the 
impact of lensing by the cluster mass. This consists of two 
effects: a magnification effect whereby the background pop- 
ulation is brightened, and a tangential compression of the 
the background volume. While the former can increase the 
counts, the latter can decrease it (i.e., the effect is to nar- 
row and lengthen the observable volume behind t he cluster) . 
The p recis e details of this correcti on are given bv lTrenthaml 
jl99S|) and lBernstein et alJ dl995Tl . 

In Fig. |H|we plot the fractional change in the observed 
number of 'field' galaxy counts caused by lensing, fi ens , as 
a function of l/-band magnitude for each cluster - averaged 
over the central 500 kpc. That is, the ratio of the number of 
'field' galaxies observed at a given magnitude viewed along 
the cluster line-of-sight to the number of galaxies expected if 
not observed through the cluster lens. In modelling the lens- 
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Figure 8. The function fi ena averaged over the central 500 kpc 
for A2218 (solid line), A2443 (dashed line) and A119 (dotted line). 

ing ef fect we follow the procedure outlined in IPracv et alJ 
(2004) and references therein. We use a value for the ve- 
locity dispersion of 1370 kms" 1 for A2218 jLe B orgne et all 
1992) and 778 kms" 1 for A119 iStruble fc Roodll987l) : since 
no velocity data is available for A2443 we assume a velocity 
dispersion for this cluster of 1000 kms -1 . We use the local 
lumin osity function measured from the MGC iDriver et alJ 
120051) and convert to the V-band using B — V = 0.94 
jNorberg et"aT1l2002l) . 

In Fig. [S]the function fi ens is shown for A2218 as the 
solid line, for A2443 as the dashed line and for A119 as 
the dotted line. Clearly, the effect is largest for A2218, pri- 
marily as a result of the cluster's higher redshift. In A2218 
the effect is most evident for galaxies with V fs 21 with 
a fractional change in the 'background' population of ap- 
proximately 10%. At this level the lensing effect represents 
a minor component of the error budget. However it is worth 
noting that the lensing correction depends critically on the 
assumptions concerning the evolution of the galaxy number 
counts as a function of redshift, which is not well known. 
Given this uncertainty and the small size of the corrections 
in Fig. |HJ we do not correct for lensing at this time. 



4 LUMINOSITY FUNCTIONS 

We now use our data along with our measured 'field' 
counts to study the luminosity distribution of galaxies in 
each cluster. The faint-end limits for A2218, A2443 and 
A119 are, respectively, My = -16.7, My = -15.4 and 
My = —13.3 mag. T hese include correct ion for galactic d ust 
jSchleeel et al.lll998l) and K-correction iPoggiantilll997l) . 

4.1 Central luminosity distribution 

We first construct, for each cluster, the luminosity distri- 
bution of galaxies within a cluster-centric radius of 1 Mpc. 
We use equation JHJ - scaled to the appropriate area - to 
correct 'statistically' for contamination by the superimposed 
'field' galaxy population. These are shown in the left panels 
of Fig. H for A2218 (top), A2443 (middle) and A119 (bot- 
tom) - and are also tabulated in Table [5] We fit each of 
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Figure 9. Left panels: The luminosity distributions for the clus- 
ters A2218 (top), A2443 (middle) and A119 (bottom) within a 
cluster-centric radius of 1 Mpc. The da ta are shown by the stars 
and the solid line is a ISchechterl <197rj) function fit to the data. 
The right panels show the 1, 2 and 3<r error ellipses for the 
Schechter function parameters M* and a. Also shown (open cir- 
cle in top panel) is the best— fitting Schechter function parameters 
for the 'bright-end ' of the luminosity distribution of A2218 from 
IPracv et aLM2004l) . 

the luminosity distributions with a ISchechterl il976Tl func- 
tion; which has been convolved with the magnitude bin 
width. The luminosity distributions in all three clusters are 
well described by a Schechter function with shape param- 
eters M v = -21.541°;^ and a = -l.U±° ° 8 7 for A2218, 
My = — 21.32l°;|g and a = — l.ll±° g for A2443 and 
My = -21.51± ;^ and a = -1.22±°;°g for A119. Since the 
parameters M* and a are dependent we plot, in the right 
hand panels of Fig. [51 the 1, 2 and 3cr error contours for the 
Schechter function fit. 

The lumino sity distribu t ion o f galaxies in A2218 has 
been studied bv IPracv et al] J20041) using a deep HST mo- 
saic of the cluster. Their 'bright-end' luminosity distribu- 
tion is comparable in bandpass, areal coverage and luminos- 
ity range to the one derived above - with a projected area 
covered of ~ 2.6 /i^ 2 Mpc 2 (compared with ~ 3.1 2 Mpc 2 
in this study) and a luminosity range corresponding to 
My < —17 (compared with My < —16.7 here). Here 
we have conver ted to the V-band using the relation V = 
F606W + 0.12 jNorberg et al.ll200llDriver et al.ll2003D and 
corrected for differences in assumed cosmology, dust ex- 
tinction, K-correction, and the omission of the cD galaxy 



Luminosity segregation in three clusters of galaxies 9 



Table 2. luminosity distributions 



Cluster 


Mag interval 


Counts 


Error 


A2218 


-23.68 


< 


M v 


< 


-22.68 


2.19 


1.73 




-22.68 


< 


My 


< 


-21.68 


26.67 


5.39 


„ 


—21.68 


< 


M v 


< 


—20.68 


60.48 


8.19 


» 


-20.68 


< 


M v 


< 


-19.68 


82.16 


10.00 




-19.68 


< 


M v 


< 


-18.68 


103.37 


12.29 




-18.68 


< 


M v 


< 


-17.68 


138.87 


16.22 


„ 


-17.68 


< 


M v 


< 


-16.68 


212.23 


22 98 


A2443 


-23.39 


< 


M v 


< 


-22.39 


0.42 


1.00 




-22.39 


< 


M v 


< 


-21.39 


16.28 


4.24 


„ 


-21.39 


< 


M v 


< 


-20.39 


31 01 


6.00 


„ 


-20.39 


< 


M v 


< 


-19.39 


47 78 


7.87 


„ 


-19.39 


< 


M v 


< 


-18.39 


66.31 


10.30 


„ 


-18.39 


< 


M v 


< 


-17.39 


100.41 


14.46 




-17.39 


< 


M v 


< 


-16.39 


99.91 


19.78 




-16.39 


< 


M v 


< 


-15.39 


73.38 


28.96 


A119 


-23.31 


< 


M v 


< 


-22.31 


1.80 


1.41 




-22.31 


< 


M v 


< 


-21.31 


12.24 


3.61 




-21.31 


< 


M v 


< 


-20.31 


32.24 


5.92 




-20.31 


< 


M v 


< 


-19.31 


66.54 


8.72 




-19.31 


< 


M v 


< 


-18.31 


63.48 


9.72 




-18.31 


< 


M v 


< 


-17.31 


65.33 


12.72 




-17.31 


< 


M v 


< 


-16.31 


111.00 


20.16 




-16.31 


< 


M v 


< 


-15.31 


173.38 


32.21 




-15.31 


< 


M v 


< 


-14.31 


200.41 


49.05 




-14.31 


< 


M v 


< 


-13.31 


254.88 


78.41 



in IPracv et alJ l)2004|l . We show the be s t-fitti ng Schechter 
function parameters from IPracv et al.l <l2004t) along with 
the error contours in Fig. [5] The recovered Schechter func- 
tion parameters are in good agreement ( — 1.23 ± 0.09 cf. 
a = — 1.14±g in this study and M v = -21.85 ±0.46 com- 
pared with M v = -21.54j£u here- 
in Table [3] we tabulate the results of recent wide-field 
V-band LF studies: column (1) identifies the cluster, column 
(2) gives the cluster redshift, column (3) gives the total area 
covered, column (4) indicates the absolute magnitude lim- 
its of the data, column (5) presents the characteristic mag- 
nitude, column (6) gives the value of the faint-end-slope 
parameter, and column (7) identifies the reference. All val- 
ues have been converted to an Qm = 0.3, = 0.7 and 
Ho = 70kms _1 Mpc _1 cosmology. Although, comparison 
of the LF parameters is somewhat complicated by differ- 
ences in redshift (and hence rest-frame bandpass), areal cov- 
erage, and the magnitude limits of the data - the values, 
and the precision, of the LF parameters presented here are 
comparable to those given in the literature. However, we 
derive faint -end-slope parameters for these clusters which 
ar e marginally flatter than the ave rage. With the exception 
of lAndreon fc Cuillandre's f2002ll Coma LF, our Abell 119 
data represents the deepest determination of a wide-field 
V-band LF available. 




Figure 10. Left panel: The composite luminosity function of the 
three clusters. Right panel: the 1, 2 and 3<r error ellipses with 
a set of Schechter function parameters for composite luminosity 
functions from the literature superimposed. See text for details. 



subtracted number counts in each magnitude bin normalised 
by the cluster's 'richness': 



C(M) = 



(10) 



where C(M) is the combined galaxy number 'counts' in 
the magnitude interval centred on M, n c is the number of 
clusters in the composite luminosity distribution (n c = 3) 
and Ni(M) is the galaxy number counts in the magnitude 
bin centred on M for the ith cluster. Ri is the 'richness 
count' of the ith cluster, which is taken to be the total 
number of background-subtracted galaxies in the magnitude 
range —21 < My < — 19 within a cluster-centric radius of 
1.5 ft^Mpc llValotto et al.ll200llh 

The composite luminosity distribution and best-fitting 
Schechter function (solid line) is shown in Fig. 1101 The shape 
of the best-fitting Schechter function is parameterised by 
M v = -21.82±°;^ and a = -1.22t ;£j. The 1, 2 and 3a 
error contours of the parameters M v and a for our com- 
posite luminosity distribution are shown in the right panel 
of Fig. 1101 along with the positions in the M v -ct plane 
of a subsample of the composite luminosity function pa- 
rameters from the literature (converted to the V-band). 
The open circle and star represent the field galaxy lumi- 
nosity functions derived from the MGC dDriver et alJl2005l) 
and Two Degree F ield Galaxy Redshift Survey (2dFGRS; 
IColless et al"ll200ll) . both of which exhibit faint-end slopes 
(a) that are similar to that of our composite luminosity 
function but have a brighter and fainter characteristic mag- 
nitude (A/*), respectively. The open square represents the 
best-fittin g Schechter function pa rameters for a composite 
of clusters iDe Propris et aDl2003l) drawn from the 2dFGRS 
with a mean redshift of z = 0.12. Our composite luminosity 
function has a marginally brighter characteristic magnitude 
to that of the 2dFGRS cluster composite and a slightly flat- 
ter a. The open t riangle shows the f ormal Schechter function 
parameters from lGoto et alJ i2002T) for a composite of clus- 
ters in the redshift range 0.2 < z < 0.25 from the Sloan 
Digital Sky Survey, which has a significantly brighter char- 
acteristic magnitude and flatter faint-end slope. 



4.2 Composite luminosity distribution 

We now construct a composite luminosity distribution from 
our three clusters. To do this we average the background- 



5 GEOMETRIC DEPROJECTION 

The wide field-of-view of the WFC mosaics provides cover- 
age of the entire cluster area, in the sense that the number 



10 Michael. B. Pracy et al. 

Table 3. Recent measurements of the LF parameters from wide— field V-band CCD photometry 



Cluster 


redshift 


Area (Mpc 2 ) 


Magnitude ran; 


5e 


M* 


a 


Reference 


Comaf 


0.023 


- 0.9 


-21.51 < M v < - 


11.51 


NA 


_-, 4 o+0.06 

-1.59 ±0.02 


Andrcon & Cuillandrc 


Coma 


0.023 


~ 1.1 


-21.51 < M v < - 


14.01 


NA 


Lobo et al. f 19971 


Abell 2151 


0.037 


~ 8.3 


M v < -15.0 




-21.561S;** 
-21.51^ 

91 o 9 +0.39 
Zl.oz_ 3g 


-1 29+ 09 
i 99 +0.07 
11 -0.06 
i 11+0.10 
L1L -0.09 


Sanchez- Jansscn et al. 


Abell 119 
Abell 2443 


0.044 
0.108 


~ 3.1 
~ 3.1 


My < -13.3 
M v < -15.4 




This work 
This work 


Abell 2218 


0.181 


~ 3.1 


M v < -16.7 




-21 54+ ' 14 


_1 14+O.O8 
^ -0.07 


This work 


ABCG 209 


0.209 


~ 3.4 


M v < -17.6 




-22.18 ±0.3 


-1.27 ±0.10 


Mercurio et al. f 20031 


ABCG 209 


0.209 


~ 6.9 


M v < -17.6 




-22.03 ±0.3 


-1.25 ±0.08 


Mercurio et al. ("20031 


AC 118 


0.310 


~ 5.5 


M v < -18.2 




-20.93 ±0.4 


-1.02 ±0.18 


Busarello et al. (20021 



f Based on a Schechter function fit to on-line data 



of galaxies per unit area no longer decreases with radius at 
the outermost regions covered b y the imaging. 

We use the schem e of iFabian et al.l il98ll see 
iBeiiersbergen et all2002T) to perform geometric deprojection 
of the clusters. The projected cluster galaxy population is 
separated into a set of concentric annuli within which the 
galaxy density is assumed to be constant. The deprojection 
scheme assumes the cluster galaxies are distributed spheri- 
cally symmetrically so that the projected number of galaxies 
in each annulus corresponds to a real galaxy number density 
in the corresponding spherical shell. The number density in 
each shell can then be corrected for the projection of galaxies 
from all shells of greater radius, beginning with the outer- 
most shell and iteratively moving inward. 

The projected galaxy density for each magnitude bin, 
m, in the ith annulus, Ni, m , is given by: 

Ni, m = 4>i, m F(Ri, Ri-i, Ri) + 

n 

4>j,m(F(Rj, Ri-i, Ri) — F(Rj-i, Ri-i, R4)) (11) 

j=i+i 

where the 7?i's are the radii of the annuli, 4>i, m is the 3D 
galaxy density and the function F is denned by: 



F(a,f3,~f) = -7TQ 3 



3/2 



-11-^ 



3/2 



(12) 



iBeiiersbergen et alJlifjoU) . 

We set radial bin partitions at cluster-centric radii 
of 0.3, 0.6 and 1.5 Mpc - splitting the cluster field into 
four regions (n = 4) - and perform geometric deprojec- 
tion of the data in 1 magnitude intervals. We first sub- 
tract from the number counts in each annuli the contribu- 
tion from the superimposed field galaxy population using 
equation JSJ - normalised to the area of 'detectability' for 
each annulus. We use the outermost annuli (R > 1.5 Mpc) 
to begin the deprojection, that is, set i — 3 in equa- 
tion lllll and then iteratively calculate the deprojected (3- 
Dimensional) luminosity function in each of the inner three 
annuli: 0.0 Mpc < r < 0.3 Mpc, 0.3 Mpc < r < 0.6 Mpc and 
0.6 Mpc < r < 1.5 Mpc. The errors in the 3D galaxy density 
are calculated by propagating the appropriate galaxy-count 
and field-subtraction errors through equation lllll . 

The deprojected 3D luminosity functions are shown in 
Fig.[TT]for A2218 {top), A2443 (middle) and A119 (bottom). 
The open circles, open squares and open triangles are the 
galaxy densities in the inner, middle and outer annuli, re- 



Table 4. Schechter function parameter values for the deprojected 
(3-dimensional) luminosity functions in three annuli: 0.0 Mpc < 
r < 0.3 Mpc, 0.3 Mpc < r < 0.6 Mpc and 0.6 Mpc < r < 1.5 Mpc 

Cluster Radius (Mpc) M v a 



A2218 


0-0.3 


99 -19+0.57 
zz.iz_ 57 


_0 94+ ' 19 
u ' y *-0.16 




0.3-0.6 


09 01 +0.60 


-1 30+ 016 




0.6-1.5 




-1 36+ 011 


A2443 


0-0.3 


on 7Q+0.5S 
zuwz -0.63 


-1 00+ ' 22 




0.3-0.6 


91 79 + 0.72 
Zi -' z -0.63 


-1 19+ 014 




0.6-1.5 


-20.03+Jl 


-1 01+ ' 45 

lul -0.18 


A119 


0-0.3 


-22.24+^g 


1 9C+O.I6 
izo -0.11 




0.3-0.6 


-21 5S+° :71 


-1 16+ 016 




0.6-1.5 


— 21 6fi + "' 51 

Z1.DD_ 54 


-1 90+O.I2 

1 - zo -o.io 



spectively, and the solid, dashed and dotted lines show their 
best fitting Schechter functions. The 1, 2 and 3a error con- 
tours for the Schechter function shape parameters My and 
a are shown in the right hand panels of Fig. 1111 The core re- 
gion (r < 300 kpc) of A2218 has a LF with a natter faint -end 
slope than the LFs measured at larger cluster-centric radii 
within this cluster, while the other clusters show no signif- 
icant systematic correlation with cluster-centric radius. A 
summary of the Schechter function fit parameters are given 
in Table H 



6 SUMMARY AND DISCUSSION 

We have exploited wide-field imaging of three galaxy clus- 
ters to examine the luminosity and spatial distribution of 
their member galaxies and to search for evidence of luminos- 
ity segregation. We cover a total projected area which cor- 
responds to 120.5 Mpc 2 , 50.4 Mpc 2 and 9.7 Mpc 2 in A2218, 
A2443 and A119, respectively. We find that the radial dis- 
tribution of galaxies can be well described by a King profile 
with a core radius which is essentially independent of lumi- 
nosity - the exception being the brightest galaxies in A2218 
which exhibit a more compact spatial distribution. Galaxies 
at different luminosities having the same spatial distribu- 
tion, of course, indicates that no luminosity segregation is 
present. We find, in all three clusters, a luminosity distri- 
bution of galaxies which is well described by a Schechter 
function with a flat faint-end slope (a ~ —1.2). We have 
performed a geometric deprojection of the cluster popula- 
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Figure 11. Left panels: the deprojected luminosity distributions 
of galaxies in cluster-centric annuli of 0-0.3 Mpc (open circles), 
0.3— 0.6Mpc (open squares) and 0.6-1.5 Mpc (open triangles) and 
their best fitting Schechter functions (solid, dashed and dotted 
lines, respectively). Right panels: 1, 2 and 3<r error contours for 
the Schechter function parameters My and a. A2218: top panel, 
A2443: middle panel and A119: bottom panel. 

tion. The luminosity function remains flat after deprojection 
and (except for the core region of A2218) we find no statisti- 
cally significant evidence for any change with cluster-centric 
radius - consistent with our radial profile analysis. 

The segregation of galaxies of different luminosities 
is a key prediction of hierarchical models of galaxy for- 
mation and evo l ution . The numerical CDM models of 
iKauffmann et all (119971) . for example, predict that low lumi- 
nosity galaxies s hould be less clustere d than the bright 'gi- 
ant' galaxies (sec Phillipos et al. 1998). Evolutionary mecha- 
nisms operating in clusters which should give rise to luminos- 
ity segregation have also been suggested. iMoore et al.l(|l998l) 
have been able to explain the dwarf population-density re- 
lation (in which the ratio of dwarf galaxies to giant galaxies 
increases with decreasing local galaxy density) as the result 
of galaxy 'harassment' - which operates more effectively in 
regions of high galaxy density. Interaction with the cluster 
tidal field will destroy the lowest surface brightness (least 
bound) objects in the central region of a cluster, resulting 
in a deficiency of faint galaxies in the cluster core. If lumi- 
nosity segregation in clusters is not present then the impact 
of these effects on the low-luminosity cluster galaxies has 
been exaggerated. 

The key issue in searching for luminosity segregation 



in clusters is that of the technique of background subtrac- 
tion employed. The fractional difference between the total 
(cluster+'field') galaxy counts and the 'field' galaxy counts 
decreases with decreasing luminosity. Therefore, any system- 
atic error in the estimate of the field galaxy counts can mimic 
both a faint-end upturn in the luminosity function and lu- 
minosity segregation - in the sense that underestimating the 
field counts will result in a larger fraction of faint galaxies 
at low local galaxy densities. Using spectroscopy to confirm 
cluster membership is essentially useless for all but the near- 
est clusters, since the galaxies at the critical faint end of the 
luminosity function are too faint for spectroscopic confirma- 
tion. Hence, most often a statistical subtraction of the field 
is performed using galaxy counts measured from a regi on off- 
set from the cluster line-of-sight. IValotto et aD l|2QQll) point 
out one possible source of bias in statistical field subtrac- 
tion due to the alignments between filaments and clusters 
which can result in an underestimate of the 'field' galaxy 
population. Here we have used an analytic representation 
(equation |SJ of the radial distribution of galaxies toward 
the cluster sight-line to simultaneously obtain the cluster 
profile and the superimposed 'field' galaxy surface density. 
Note, fitting equation © to the radial profiles in Fig. [3] 
Fig. 0] and Fig. does not require prior subtraction of the 
field galaxy population and we, therefore, suggest that this 
is a superior test for luminosity segregation. Very recently 
lAndreon et alJ l)20flftl have suggested an alternative method 
to alleviate and quantify the problems of background sub- 
traction in deriving cluster luminosity functions - a detailed 
comparison of these methods and the traditional method 
should be pursued in a future paper. 
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